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Abstract 

We investigate the origin of the cos 2(/> azimuthal asymmetry in unpolarized semi- 
inclusive DIS. The contributions to this asymmetry arising from the intrinsic trans¬ 
verse motion of quarks are explicitly evaluated, and predictions for the HERMES 
and COMPASS kinematic regimes are presented. We show that the effect of the 
leading-twist Boer-Mulders function hj^{x,\i.‘^), which describes a correlation be¬ 
tween the transverse momentum and the transverse spin of quarks, is quite signifi¬ 
cant and may also account for a part of the cos 2(f> asymmetry measured by ZEUS 
in the perturbative domain. 
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1 Introduction 


The importance of the transverse-momentum distributions of quarks for a full 
understanding of the structure of hadrons has been widely recognized in the 
last decade [1,2,3,4]. In semi-inclusive deep inelastic scattering (SIDIS), the 
k^-dependent distributions give rise to various azimuthal and/or single-spin 
asymmetries, which are currently under direct experimental scrutiny [5,6]. 
Two leading-twist distributions of great relevance for their phenomenological 
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implications are the Sivers function f^{x, kf,) [7] and its chirally-odd part¬ 
ner the so-called Boer-Mulders function [4], These two distribu¬ 

tions describe time-reversal odd correlations between the intrinsic momenta 
of quarks and transverse spin vectors [8]. In particular, represents an 
azimuthal asymmetry of unpolarized quarks inside a transversely polarized 
hadron, whereas represents a transverse-polarization asymmetry of quarks 
inside an unpolarized hadron. Recently, it has been proven by a direct cal¬ 
culation [9] that and hf are non-vanishing: interference diagrams with a 
gluon exchanged between the struck quark and the target remnant generate 
non-zero asymmetries. The presence of a quark transverse momentum smaller 
than Q ensures that these asymmetries are proportional to M/Zc^, rather than 
to M/Q, and therefore are leading-twist quantities. Moreover, a careful con¬ 
sideration of the Wilson-line structure of k^-dependent parton densities shows 
that /fi^ and are not forbidden by time-reversal invariance [10,11] (for a 
possible chiral origin of these distributions, see [12]). 

The Sivers function is known to be responsible for a sin((;/) — (ps) single¬ 
spin asymmetry in transversely polarized SIDIS [5,6,13]. The Boer-Mulders 
function hf produces azimuthal asymmetries in unpolarized reactions. Boer 
[14] argued that it can account for the observed cos 20 asymmetries in unpo¬ 
larized 71N Drell-Yan processes [15,16]. This was quantitatively conhrmed in 
[17,18], where hf was calculated in a simple quark-spectator model and shown 
to explain the Drell-Yan data fairly well. 

A similar cos 20 asymmetry occurs in unpolarized leptoproduction. As we 
shall see, there are three possible mechanisms generating this asymmetry: 
1) non-collinear kinematics at order k^/Q"^ [19]; 2) the leading-twist Boer- 
Mulders function [4] coupling to a specular fragmentation function, the so- 
called Collins function [20], which describes the fragmentation of transversely 
polarized quarks into unpolarized hadrons; 3) perturbative gluon radiation 
[21,22,23,24]. The purpose of this paper is to study the first two sources of the 
cos 20 asymmetry, both related to the intrinsic transverse motion of quarks. 
They are especially relevant in the HERMES kinematic regime ((Q^) ~ 2 
GeV^), but the Boer-Mulders contribution, being leading twist, can also sur¬ 
vive at higher and partly account for the asymmetry measured by ZEUS 
in this domain [25]. 

In recent years, the cos 20 asymmetry in leptoproduction was phenomenologi¬ 
cally studied by some authors [26,27]. In [26] only the 0{k‘^/Q‘^) term and the 
perturbative contribution were included, whereas the Boer-Mulders effect was 
not considered. Our calculation is more similar to that presented in [27], the 
main differences being that we use a model for adjusted on the Drell-Yan 
data [18], and compute the asymmetry according to its experimental dehnition 
(which incorporates a cutoff on the transverse momentum of the hnal hadron). 
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Fig. 1. Lepton and hadron planes in semi-inclusive deep inelastic scattering. 

2 The cos 20 asymmetry in unpolarized SIDIS 


The process we are interested in is unpolarized SIDIS: 


m + p(p) ^ i\t) + h{Pn) + x{Px) . 


(1) 


The SIDIS cross section is expressed in terms of the invariants 

g' P-q P-Ph 

2P ■ q P ■ t P-q 


( 2 ) 


where q = £ — £' and We adopt a reference frame such that the 

virtual photon and the target proton are collinear and directed along the z 
axis, with the photon moving in the positive z direction (Fig. 1). We denote 
by Ict’ the transverse momentum of the quark inside the proton, and by 
the transverse momentum of the hadron h. The transverse momentum of h 
with respect to the direction of the fragmenting quark will be called p^. All 
azimuthal angles are referred to the lepton scattering plane (we call 0 the 
azimuthal angle of the hadron h, see Fig. 1). 

Taking the intrinsic motion of quarl^s into account, the SIDIS cross section 
reads at leading order 


da 

dxdg dz d^Pr 


27raLg 


Y,elx[l + {1 - yf] 


X y d^kr J - zkr - Pt) D^{z,Pt) ^ 


(3) 


where fi{x, k|) is the unintegrated number density of quarks of flavor a and 
PH^iPt) is the transverse-momentum dependent fragmentation function of 
quark a into the hnal hadron. We recall that the non-collinear factorization 
theorem for SIDIS has been recently proven by Ji, Ma and Yuan [28] for 

Ft <g. 
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As shown long time ago by Cahn [19], the transverse-momentum kinematics 
generates a cos 20 contribution to the unpolarized SIDIS cross section, which 
has the form 


X J d^kr J d^pT 5 ^(Pt — 2:kT — Pt) 

X - — fi{x,\4) D^{z,p^t) cos 20, (4) 

where h = Pt/Pt- Notice that this contribution is of order k^/Q'^, hence it 
is a (kinematic) higher twist effect. 

The second dependent source of the cos 20 asymmetry involves the Boer- 
Mulders distribution h'l; coupled to the Collins fragmentation function of 
the produced hadron. The explicit expression of this contribution to the cross 
section is [4] 


daPP 


dx dy dz d^Pr 




J2elx{l-y) 


COS 24 > ^ ^ 

2l, f X2/ 


X y d^kr y d^pT 5 ^(Pt — 2:k'r — Pt) 

2 h ■ kr h ■ pt - kr ■ Pt ,_±a/ 


X 


h0“(x, k^) p^) cos20 , (5) 


zMMh 


It should be noticed that this is a leading-twist contribution, not suppressed 
by inverse powers of Q. 

The asymmetry measured in experiments is dehned as 


(cos 20) = 


/ da cos 20 
J da 


(6) 


where the integrations are performed over the measured ranges of x, y, z and 
with a lower cutoff Pc on Pt, which is the minimum value of Pt of the detected 
charged particles. Using Eqs. (3) and (5), (cos20/i) is given by 


(cos 20) 


Jiff Eg e! 2x(l - y) {A[f^, D<^] + ^ 

I!IIEaelx[l + {l-yy]C[n,Df] 


(7) 
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where 


^T',max 


X2 J/2 22 


Illh / dP.P./dx/d,/d. 


XI yi zi 


and (x is the angle between and k^) 


( 8 ) 


^[/“, D^] = J d^kr J d^pT 5^(Pr - zkr - pr) 

2 (kr ■ h)^ — k^ 


X 


ra 


2it 


= / dkT kx / dx 


/f(a;,k^)L)“(x,Pr) cos 20 
2k|, cos^ X ~ 




x/f(a;,k|)D“(x,|PT-;^kTr), 


(9) 


= j d^kr j d^PT S\Pt - zkT - Pt) 

2 h ■ kr h ■ pt - kr ■ pr ,±a 


X 


^r(a^,kT)^i “(^,Pr) 


zMMh 

= /dx‘‘t + (PrX)<^TCosx-2k|cos^x 


0 0 

X k^) |Pr - ;2k 


MMh 

| 2 \ 


( 10 ) 


c[/f,D;|=|d=kT |dXT^"(PT-^k7.-PT)/f(I,k?,)D;(^,p?,) 

(50 27r 

= J dkrkr J dx/r(x,k^)D“(x,|Pr-;^krP). (11) 

0 0 


3 Calculation and Results 


In order to calculate (cos 20) one needs to know the kx- and pT-dependent dis¬ 
tribution and fragmentation functions appearing in Eqs. (9)-(ll). Independent 
information on the Boer-Mulders function /i0(a;, k|,) can be obtained from the 
study of the cos 20 azimuthal asymmetry in unpolarized Drell-Yan processes, 
which has been measured in ttN collisions [15,16]. In [17,18] this asymmetry 
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Fig. 2. The SIDIS cos 2(j) azimuthal asymmetry in the HERMES domain as a func¬ 
tion of the cutoff Pc, for two values of The dotted curve is the leading-twist 

Boer-Mulders contribution, the dashed curve is the higher-twist term, the solid 
curve is the sum of the two contributions. 

was estimated by computing the distribution of the pion and of the nu¬ 
cleon in a quark spectator model [29,30]. To compute the cos 20 azimuthal 
asymmetry in SIDIS we adopt the same distributions h^i^x, kf.) and /i(x, k^) 
used in [18]. We assume that the observables are dominated by u quarks (i.e., 
we consider vr’*' production). The set of the transverse-momentum dependent 
distribution functions is (for simplicity, we consider a spectator scalar diquark 
[30,18]) 


= N{l-xf 


{xM -I- mY + k|, 

{L^ + ’ 

,3 M {xM + m) 


( 12 ) 

(13) 


where is a normalization constant, m is the constituent quark mass, and 
L? = {1 — x) + X — x{l — x) . (14) 


Here A is a cutoff appearing in the nucleon-quark-diquark vertex and Md is the 
mass of the scalar diquark. As it is typical of all model calculations of quark 
distribution functions, we expect that Eqs. (12) and (13) should be valid at low 
values, of order of 1 GeV^. The average transverse momentum of quarks 
inside the target, as computed from ( 12 ), turns out to be — 0.54 GeV. 

Goming to the fragmentation functions, for we adopt the simple parametriza- 
tion suggested by Gollins [20] 

McMn 
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Fig. 3. The SIDIS cos 2(j) azimuthal asymmetry in the HERMES domain, as a func¬ 
tion of X with Pc = 0.5 GeV. The dotted curve is the leading-twist Boer-Mulders 
contribution, the dashed curve is the higher-twist term, the solid curve is the sum 
of the two contributions. 

where Mq is a free parameter. We assume a Gaussian dependence for the 
unintegrated unpolarized fragmentation function: 

D,(^,p|) = r>iW^e-'>W), ( 16 ) 

7r{p^) 


so that / d^pT Di{z, pf.) = Di{z). Finally, the integrated unpolarized fragmen¬ 
tation function for pions Di{z) is taken from the Kretzer-Leader-Christova 
parametrization [31], 

Di{z) = 0.689 (1 - (17) 


valid at = 2.5 GeV^. For the parameters in Eqs. (12) and (13) we choose 
the values Md = 0.8 GeV, m = 0.3 GeV, A = 0.6 GeV, as = 0.3, which 
are the same as in [18]. As for the parameters in Eqs. (15) and (16), we £x 
Me to 0.3 GeV and show results for two values of the average transverse 
momentum: = 0.5 GeV and 0.6 GeV (we checked that a variation of 

Me is reproduced by a change of {pY)^^Y- 

The HERMES kinematics is characterized by the following ranges: 0.02 < 
X < 0.4, 0.1 < y < 0.85, 0.2 < 2 ; < 1, {QY = 2 GeV^. Our predictions 
for the cos 20 asymmetry in this regime are displayed in Fig. 2, where we 
show separately the higher-twist term and the leading-twist Boer-Mulders 
contribution. For a typical transverse momentum cutoff Pc = 0.5 GeV, these 
two terms are comparable and the predicted asymmetry lies in the range 
(cos 20) = 0.02—0.04. The x-dependence (with z integrated over the accessible 
interval) and the z-dependence (with x integrated over the accessible interval) 
are shown in Figs. 3 and 4, respectively. As one can see, the asymmetry is 
larger at small x and large z. 
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Fig. 4. The SIDIS cos 2(j) azimuthal asymmetry in the HERMES domain, as a func¬ 
tion of z with Pc = 0.5 GeV. The dotted curve is the leading-twist Boer-Mulders 
contribution, the dashed curve is the higher-twist term, the solid curve is the sum 
of the two contributions. 



Eig. 5. The SIDIS cos 2(j) azimuthal asymmetry in the COMPASS domain, as a 
function of x with Pc = 0.5 GeV. The dotted curve is the leading-twist Boer-Mulders 
contribution, the dashed curve is the higher-twist term, the solid curve is the sum 
of the two contributions. 

In Fig. 5 we plot our results for the x-dependent asymmetry (integrated over 
z) in the COMPASS kinematic domain. The correlation between x and is 
such that the lowest x bin {x = 0.005) corresponds to ~ 1 GeV^, whereas 
the highest x bin in Fig. 5 (a; = 0.25) corresponds to ~ 24 GeV^. Again, 

the asymmetry is of order of few percent and decreases with x. 

There are available data on the cos 20 asymmetry in SIDIS coming from the 
ZEUS experiment [25]. The ZEUS kinematic ranges are: 0.01 < x < 0.1, 
0.2 < y < 0.8, 0.2 < 2 : < 1, > 180 GeV^. At such large values, the 















Pc GeV 

Fig. 6. The SIDIS cos 2(f> azimuthal asymmetry as a function of the cutoff Pc in the 
ZEUS domain. Data are from [25]. 


higher twist contribution is clearly irrelevant. Since only the evolution of 
the kT moments of is known [32], and not that of itself, we assume for 
simplicity that the distributions (12) and (13) scale exactly, i.e. that they are 
valid for any (one should recall however that Sudakov form factors aris¬ 
ing from soft gluon contributions may reduce the Boer-Mulders asymmetry at 
very high [33]). The result for the cos 2(f) asymmetry in the ZEUS kinematic 
domain is shown in Fig. 6, where it is compared with the experimental data. 
The agreement is rather good for low values of the Pt cutoff (up to 0.5 GeV). 
For larger Pt values one expects of course a relevant perturbative contribu¬ 
tion. Including this contribution is beyond the purpose of this paper, which 
is primarily devoted to predictions for the low-Q^ domain. A more extended 
analysis of the cos 20 asymmetries, taking into account also the perturbative 
term, is in progress and will be reported soon [34]. 

For completeness we recall that long time ago the European Muon Collab¬ 
oration at CERN measured (cos 20) for > 4 GeV^ [35]. The EMC data, 
however, are affected by large uncertainties and do not allow drawing deh- 
nite conclusions about the magnitude and the shape of the asymmetry. The 
comparison of our predictions with these data is shown in Fig. 7. 

In conclusion, we predicted the cos 20 asymmetry for semi-inclusive deep in¬ 
elastic scattering in the kinematic regions of the HERMES and COMPASS 
experiments. We found that (cos 20) is of order of few percent and tends to 
be larger in the small-x and large- 2 ; region. The combined analysis of the fu¬ 
ture data on (cos 20) and of the previous ZEUS measurements in the high-Q^ 
domain (where higher twist effects are irrelevant) will allow to get informa¬ 
tion on the Boer-Mulders function, shedding light on the correlations between 
transverse spin and transverse momenta of quarks. 
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Fig. 7. The cos 20 azimuthal asymmetry (divided by/ 2 (y) = (1— 2 /)/[l+(l—y)^]) as a 
function of xp = 2Pp/W as measured by EMC [35]. The curves are our predictions. 
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